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ABSTRACT 

We have obtained a deep, sub-arcsecond resolution X-ray image of the nuclear 
region of the luminous galaxy merger NGC 6240 with Chandra, which resolves 
the X-ray emission from the pair of active nuclei and the diffuse hot gas in 
great detail. We detect extended hard X-ray emission from kT ~ 6 keV (~70 
million K) hot gas over a spatial scale of 5 kpc, indicating the presence of fast 
shocks with velocity of ~2200 km s~^. For the first time we obtain the spatial 
distribution of this highly ionized gas emitting Fe XXV, which shows a remarkable 
correspondence to the large scale morphology of H2(l-0) S(l) line emission and 
Ha filaments. Propagation of fast shocks originated in the starburst driven wind 
into the ambient dense gas can account for this morphological correspondence. 
With an observed Lo.s-skoV = 5.3 x 10"^^ erg s~^, the diffuse hard X-ray emission 
is ~100 times more luminous than that observed in the classic starburst galaxy 
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M82. Assuming a filling factor of 1% for the 70 MK temperature gas, we estimate 
its total mass (Mhot = 1-8 x IO'^Mq) and thermal energy (Eth = 6.5 x 10^'' 
ergs). The total iron mass in the highly ionized plasma is Mpe = 4.6 x lO^M©. 
Both the energetics and the iron mass in the hot g consistent with the 

expected injection from the supernovae explosion during the starburst that is 
commensurate with its high star formation rate. No evidence for fiuorescent Fe 
I emission is found in the CO filament connecting the two nuclei. 

Subject headings: X-rays: galaxies — galaxies: starburst — galaxies: interactions 
— galaxies: ISM — galaxies: individual (NGC 6240) 



1. Introduction 



In the local universe, NGC 6240 (^=0.02448 ± 0.00003; IPownes et al.lll993l ) is a unique 
galaxy, in the throes o f a violent raerging event and on its way to becoming an elliptical galaxy 



(iTacconi et al.l Il999l : iTecza et al.l I2OOOI : iBush et all I2OO8I : Engel et all I2OIOI : Medhng et al. 



2OIII ). It is experiencing intense star formation (e.g., 61 iSOM© yr ^ in Yun & Carilli 2002; 



25M0 yr-i in Engel et al. 2010). With Lpm ~ lO^^-^L© just below IO^^^q, the luminosity 
threshold for the ultraluminous infrared galaxies (ULIRGs, Sanders & Mirabel 1996; Genzel 
et al. 1998), NGC 6240 is e xpected to becora e an ULIRG when the galaxies coalesce and a 
final starburst is triggered (lEngel et al.ll2010l ). 



The Chandra X-ray images (ILira et al.l l2002l : iKomossa et al.l |2003| ) of the central re- 



gion of NGC 6240 revealed two gravitationally interacting active supermassive black holes 

(SMBHs) separated by 1.5" (~0.7 kpc). These sources are characterized by th e highly ab- 

sorbed hard X-ray spec tra typical of Compton thick active galactic nuclei (AGNs; llwasawa fc Comastri 



19981 : iMatt et al.ll2000l ). They each have an observed luminosity of L2-8keV ~ 10^^ erg s 



and show prominent neutral Fe Ka lines at 6.4 keV present in the spectra of both nuclei 
(IKomossa et al.ll2003l ). 



Limited by the number of counts, marginal evidence of an H-like Fe line i n the spectral fit 



and ex tended hard X-ray emission in the morphology was also suggested by IKomossa et al. 



(I2OO3I ). This line emission could be associated with a spatially resolved, multi-temperature 
hot gas outfiow, powered by the sta rburst activi t y tha t dominates the 0.5-3 keV X-ray 
emission (Ln.i-2kRV 



2003 



Huo et al 



]^q43 g 1. 



Schulz et al.1 Il998l : IKomossa et"aD IgOOgL iPtak et al.l 



2OO4I : ICrimes et alllioosl ). XMM- Newton observations jBoUer et al.l I2OO3I : 



Netzer et al.l l2005l ) resolved the Fe K line complex into three narrow lines, t he neutral Fe 



Ka, Fe XXV, and a blend of Fe XXVI with Fe K/3. Using a three-zone model, iNetzer et al. 
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( 120051 ) further identified that the emission hnes of higher ionization ions (e.g., Fe XXV hne) 
originated from the inner 2.1 kpc region. These XMM- Newton data, however, do not have 
the resolution to resolve spatially the line-emitting regions. 

In this paper, we focus on the spatially resolved 5.5-8 keV X-ray emission of the in- 
nermost central ~5 kpc (~ 10") region in NGC 6240 based on new, and archival, Chandra 
imaging. This choice of hard energy range allows us to investigate the origin of the high 
temperature gas and the highly ionized iron emission. A detailed study of t he extended soft 



X-ray halo has been presented in a companion paper (INardini et al.l 120131 ). A third paper 
(Wang et al. in preparation) will discuss the extended soft X-ray emission in the r < 15" 
region. Taking advantage of the order of magnitude deeper imaging, along with subpixel 
resolution, we are able to firmly establish the presence and the spatial distribution of the 
hard X-ray emission. Throughout this work we adopt a luminosity distance of Dl = 107 
Mpc to NGC 6240 (1" = 492 pc), based on the concordance cosmological parameters {Hq = 
70.5 km s^i Mpc-i, ^]„=0.27, Qa=0-7^; Komatsu et al. 2011). 



Observations and Data Reduction 



NGC 6240 was observed on May 31, 2011 with the Chandra Advanced CCD Imaging 
Spectrometer (ACIS; ICarmire et aDl2003h for 146 ks (ObsID 12713, PI: Fabbiano). The bi- 
nary nucleus of NGC 6240 was placed near the nominal aimpoint of the backside-illuminated, 
low energy sensitive S3 chip. The ACIS data were reprocessecE] with the Chandra Interactive 
Analysis of Observations software package (CIAO; version 4.4) using the latest Chandra Cal- 
ibration Database (CALDB) version 4.5.1. The net exposure time was 144 ks after screening 
for brief periods of ele vated background. T he archival ACIS S3 observation (ObsID 1590, 
previously presented in lKomossa et al.ll2003l ) shown in this work was taken on February 22, 
2002 with a total exposure of 37 ks. We reprocessed the ACIS data to generate a new level 
2 event file th at has the latest calibration and the subpixel resolution algorithm ("EDSER", 
Li et al.ll2003l ) applied. Because the effective area of the Chandra High Energy Transmission 
Grating (HETG) zeroth-order image is about 60% of ACIS-S at 6.4 keV, we also make use 
of two archival deep HETG observations (PI: Canizares; ObsID 6908, 159 ks, observed on 
May 16, 2006; ObsID 6909, 143 ks, observed on May 11, 2006), providing an equivalent of 
180 ks ACIS-S imaging in the Fe K band. Therefore, the combined depth of our data reaches 
~363 ks, which improves by an order of magnitude over the first ACIS image of NGC 6240 
(Komossa et al. 2003) in the 5.5-8 keV range. Table [1] summarizes the Chandra ACIS and 



^ ,http : //cxc . harvard . edu/ciao/ threads/createL2/ 
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HETG observations used in our work. 



3. Imaging Spectroscopic Analysis of the NGC 6240 Central Region 

Figured] shows the Chandra ACIS image of the inner 25" x25" (~12 kpc across) region of 
NGC 6240 in the 5.5-8 keV band extracted from the merged data, together with an archival 
HST ACS /WFC F814W image (Cycle 14 program 10592, PI: A. Evans), highlighting the 
highly disturbed optical morphology of NGC 6240. The two nuclei are resolved as the 
brightest X-ray point sources in the center, with peak emission separated by 1.5" (710 pc), 
matching the nuclear compact radio source positions (e.g., Nl and N2 in the 8.4 GHz map 
of Colbert et al. 1994; Nl and S in the 1.4 and 5 GHz maps of Beswick et al. 2001; North 
and South black hole positions in Max et al. 2007). 

Diffuse emission is clearly present (Figure [T]d), extending to r = 10" (~5 kpc). Hereafter 
we refer this region as the central region of NGC 6240. The X-ray morphology resembles a 
"heart" with loop-like features in the south, the north-east, and the north-west directions 
(Figure Wp)- A comparison between the observed surface brightness profile and that of 
a point source further confirms the presence of significant hard X-ray emission reaching 
r ~ 10" radial distance. Figure |2] shows the excess of hard X-rays over the scattered emission 
in the wings of the point spread function (PSF). The radial profile flattens past 10", where 
background emission dominates. The Chandra PSF shown in Figure|2]was created simulating 
a point source at the same position on the detector as the center of NGC 6240 in the single 
observation of the longest exposure (ObsID 12713) using the Chandra Ray Tracer (ChaRTn) 
and the MARX§| simulator. This PSF was then rescaled to match the combined 5.5-8 keV 
counts from the two nuclei. Since the nuclei are separated by 1.5", and we are interested in 
demonstrating the presence of extended emission several arcseconds from the nuclei, this is 
an adequate procedure. For the investigation of extended emission close to the nuclei (e.g., 
between the nuclei), we have performed a more complex 2-D image fitting with two PSF 
components, described in the next section. 



' http : //cxc . haxvard . edu/ chart/ 



^Version 4.4.0; 



http : / / space . mit . edu/ cxc/marx/ 



- 5 - 



3.1. Spectral Analysis: Neutral, He-like and H-like Fe K lines 

Previous work by Boiler et al. (2003) identified ionized Fe K emission lines from the 
presence of a coUisionally ionized plasma with a temperature kT = 5.5 ± 1.5 keV, producing 
the Fe XXV (at 6.7 keV) and Fe XXVI (at 6.93 keV) Ka recombination lines. Boiler et al. 
also inserted in their fit the neutral or low ionization Fe Ka line (expected at 6.4 keV) and 
K/3 line (at 7.058 keV), each modelled with a Gaussian line. Such fluorescent lines of large 
(up to few keV) equivalent widths (EWs) arise from large column density gas which strongly 
attenuates directly viewed X-rays from the AGNs. 

To shed light on the spectral properties of the spatially resolved central region, we 
extracted the spectra of both nuclei and all the emission in the 10"-radius circle. Figure |3^ 
compares the ACIS spectra of the central region of NGC 6240 (circle in Figured], all emission 
within r=10", ~5 kpc from the nuclei) and the innermost region containing the two nuclei 
(ellipse in Figure [H full lengths of the axes 4" x 2"). 



The complex Fe K line structure detected in AMM- Newton spectra (IBoller et al.l 12003 



Netzer et al.ll2005l ) is clearly present in the Chandra ACIS spectra (Figure [3]) and contributes 
most counts, in contrast to the low continuum in the same region. Approximately 80% of 
this continuum comes from the innermost region enclosing the nuclei. There are notable 
differences between the two spectra, which are f urther illustrated i n Figure [3] showing the 



residual spectrum excluding the nuclei. Following iBoUer et al.l ( 120031 ). we added four narrow 
Gaussian emission lines (representing neutral Fe Ka, Fe XXV, Fe XXVI, and Fe K/3) to 
a highly absorbed power law continuum to fit the spatially resolved spectra (the northern 
nucleus, the southern nucleus, and the extended emission within r = 10", i.e., nuclei+diffuse). 
This phenomenological approach is adopted to identify the spatial distribution of the line 
emission, and a more physically meaningful model will be used in the later section. To 
prevent oversubtraction of diffuse emission from the hot gas in the nuclear region, we did 
not exclude the innermost region when evaluating the hot gas content in the r < 10" emission; 
instead, we included nuclear components to represent the contribution from the AGNs. The 
Fe XXVI and K/3 lines are blended together at the spectral resolution of the ACIS CCD, 
therefore we fixed the K/3 line center energy at 7.058 k eV (rest frame) and scaled the line 



intensity according to the theoretical ratio (1/8.8, e.g., iPalmeri et al.ll2003[ ) relative to the 
Ka line. The best-fit centroid energies, EWs, and line fiuxes, and associated 1-a errors are 
summarized in Table [2J 

The line fiuxes further confirm that the Fe XXV emission is spatially extended while 
the neutral Fe Ka emission is concentrated at the nuclei. The neutral or weakly ionized 
Fe Ka line emission at 6.4 keV (centered at 6.25 keV in the observer frame) is also mostly 
associated with the nuclei, with ~15% of the total Fe I Ka emission in the extended region 
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(Figure [3]). At the 6.4 keV energy, the PSF scattering of nuclear emission to the larger radii 
cannot be neglected, which is ~10%. Taking this into account, the spatially extended Fe 
I Ka emission is at 5% of the total. The Fe XXVI emission is fully accounted for by the 
innermost nuclear emission, indicating that the mostly ionized H-like iron is confined in this 
region. In contrast, ~30% of Fe XXV line emission appears to originate in the extended 
region outside of the nuclear region (taking into account of the PSF contribution). In the 
following we will in vestigate the morphology of this emission. Given that the XMM spectra 
( iNetzer et al.ll2005l ) come from a larger region, the similarity of the fluxes with our Chandra 
measurements excludes significant Fe emission beyond r = 10". 



3.2. Narrovi^-band Hard X-ray Images 

The analysis above shows that the AGNs alone in NGC 6240 cannot account for the 
observed Fe K complex, except for the Fe Ka and K/3 emission (the latter fainter and blended 
with Fe XXVI), which is due to fluorescence of neutral or low ionization material near the 
AGNs. Guided by the spectra (Figures [3] and H]), we extracted Chandra narrow-band line 
images to investigate the morphology of the emission line gas. 

Figure |5] shows images in the 5.5-6 keV continuum band that is free from any significant 
line emission, the Fe Ka line (6-6.4 keV) and the Fe XXV line (6.4-6.7 keV). The blended 
Fe XXVI and Fe I K/3 emission was examined but not shown here, since it does not show any 
evidence for emission beyond the nuclear region. Although the line emission dominates in 
these narrow bands, the contribution of the continuum has been removed from the images 
regardless, assuming a single power law with a photon index (F = 0.8) from the spectral fit 
of the continuum in the central region. The contribution from the continuum in the emission 
line bands is calculated using the photon index and the counts in 5.5-6.0 keV, taking into 
account the difference i n the effective area s. Given that the starburst likely dominates the 



extended emission (e.g., iNetzer et al.ll2005[ ). we further confirmed that the following results 
do not change adopting a thermal bremsstrahlung spectrum (kT~6 keV for the hot gas) to 
approximate the fiat continuum. 

The continuum-subtracted line emission images clearly show that the 6.4 keV Fe I Ka 
line peaks at the dual nuclei, with weak emission extending to r = 3" from the nuclei 
and contributing ~10% of the total ~900 Fe Ka counts. The Fe XXV emission is also 
centrally peaked at the nuclear region with a dominating southern nucleus, but it appears 
more extended to the southwest and to the east; the r = 10" region outside of the nuclei 
contains ~40% of the ~500 total Fe XXV counts. These values are in good agreement with 
the results from spectral analysis above. 
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To extract more morphological information from the data, we performed ima ge restora- 



tion using the Expectation t hrough Markov Chain Monte Carlo (EMC2) algorithm (lEsch et al. 



2004 : iKarovska et al.l l2005l ). The effectiveness of this method with Chandra imag e s was 



demonstrated with a nu mber of astronomical imaging studies ( IKarovska et al.l 12005 



2007; 



Wang et al.l l2009al 120121 ) . The EMC2 deconvolution is more ef fective for p reserving diffuse 



features compared to other PSF-deconvolution algorithms (e.g.. lLucylll974j ). Only the single 
deep imaging observation (ObsID 12713) and the associated PSF were used to minimize 
the PSF variations across multiple observations. The PSF-deconvolved images, shown in 
Figure O provide the highest resolution view of the central region. In particular, the faint 
diffuse Fe XXV emission is more obvious in Figure |6t. 



3.3. Sub-kpc-scale Extended Emission between the Nuclei 

Close examination of the images suggests the possibihty of some extended hard X-ray 
emission in the region between the two nuclei, especially in the Fe XXV band (and perhaps in 
the Fe Ka band; see Figures [5] and [6]) . However, evaluating the significance of this emission 
is very difficult because of the close separation (1.5") of the nuclei. The PSF wings of the 
two X-ray point sources overlap in this region and enhance the emission here, which needs to 
be taken into account. Simulating two point sources convolved with the PSF, we performed 
2-D image fitting to subtract the PSfQ on the 5.5-8 keV image, as well as individual narrow 
band images, using Sherpa0. The results are shown in Figure [71 For the 5.5-8 keV image 
that has most counts, besides the clear extension to the southwest and the northeast (already 
noted above as the extended central region), the residual image further indicates that there 
is extended emission in the region between the two nuclei, immediate north of the southern 
nucleus, at ~10 a significance level when evaluated using a 1" x 1" box region. For the 
narrow band images, the statistics are poorer; nevertheless, this extended feature is seen in 
all three residual images. Its significance ranges from ~6 a in the Fe XXV line emission to 
~3 a in the Fe Ka line emission. 

We further checked for possible contamination from the PSF artifaci0 found in ACIS 
data. Using the CIAO tool make_psf _asymmetry_region, this artifact should lie in position 
angle — 50° in our data. Given that this PSF asymmetry feature may contain only 5% of 



*See detailed procedure available at |http : //cxc . harvard . ed.u/sherpa4 . 4/thread.s/2dpsf / 
^Version 2 in CIAO 4.4; [http : //cxc . harvard . edu/sherpa/ 



http : //cxc . harvard . edu/ciao4 . 4/caveats/psf _cirtif act . html 
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the flujfl, it would have in any case minimal impact on our results. 



3.4. Comparison with Multiwavelength Images 

In this section, we compare the hard X-ray morphology of the center of NGC 6240 with a 
number of high spatial resolution images available in other wavelength ranges, aiming to find 
clues for the physical origin of the extended X-ray emitting gas. The absolute astrometry of 
our Chandra ACTS image is excellent, indicated by the agreement between the coordinates of 
the X-ray and the radio nuclei (~0.2" accuracy, see Max et al. 2007 for a further astrometry 
discussion). 

We first compare the 5.5-8 k eV emission wit h the optical Ha emission morphology 
from if5'r/WFPC2 F673N image Jlira et al.ll2002l ). as shown in Figure EJi. The 5.5-8 keV 
image was adaptively smoothed with CIAO tool dmimgadapt using a Gaussian kernel with 
smoothing scales of 1 pixel to 20 pixels, to better visualize structures in the diffuse emission. 
As a separate approach, the same 5.5-8 keV image with the PSF deconvolution applied is 
also shown in Figure [8b . The optic al peaks are no t align ed with the X-ray nuclei due to 
obscuration. iLira et al.l ( 120021 ) and iKomossa et al.l (|2003[ ) have previously shown that on 
10" (4.9 kpc) scale, the Ha emission and the soft X-ray emission share a similar morpholgy 
consisting of filaments and loops. Such morphological correspondence is seen but less obvious 
in the hard X-ray emission, except for the bright Ha filament to the east and the loop to the 
northwest of the nuclei. Figure [8t and [Hi further compare the Ha morphology with images in 
the 5.5-6 keV continuum band and the Fe XXV band. A close morphological correspondence 
between the Ha and the hard X-rays can be found at the base of an incomplete loop to 
the west of the nuclei (Figure [Ht and d), and at the filament to the south and east of the 
southern nucleus (Figure [8]i). 



Next we compare the 5.5-8 keV image with the [Fell] 1.644 /xm line emission (Ivan der Werf et al. 



19931 ) in Figure [Hi The [Fell] emission in NGC 6240 could be produced by shocks in supernova 
remnants (SNRs; see van der Werf et al. 1993 for detailed discussion on possible excitation 
mechanisms), and is br ightest at the loca tions of the two nuclei, consistent with the presence 
of nuclear starbursts (lEngel et al.ll2010l ). The two [Fell] emission peaks fall close to the 
X-ray nuclei, with a slight offset (0.2") to the north relative to the X-ray northern nucleus. 

We find a close correlation of the spatially extended Fe XXV emission with the near 



^See document written by Vinay Kashyap 
,http : //cxc . harvard . edu/ cal/Hrc/PSF/acis_psf _2010oct . html 



(CXC) 



available 
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IR H2(l-0) 2.12 nm emi s sion obtained w ith Keck II (IMax et al.l l2005l ). Previous work 
dvan der Werf et al.l E993I : iMax et al.l [2005! ) noted that the H2(l-0) emission in NGC 6240 
does not follow th e stellar light but peaks between the two nuclei, w hich is thought to be 
excited by shocks (IMoorwood &: 01ivalll988l : Ivan der Werf et al.lll993l ). To compare H2 and 
our X-ray images, the astrometry of the H2 image was first registered using an IR-bright 
star (2MASS J16530119+0224138, K=WA mag) in the field. Figure HO] shows a remark- 
ably good correspondence between the observed (top panel), and especially the deconvolved 
(lower panel) Fe XXV line emission and the H2(l-0) 2.12 fim emission. 

Finally we compare the Fe I Ka emission with the cold molecular gas distribution in 
the nuclear region, traced by the CO (3-2) emission obtained with the Submillimeter Array 
(SMA) from U et al. (2011). Figures [TTk and [TTb show that the CO emission resembles 
a 'bridge' or filament between the two nuclei, with the bright peak offset to the north of 
the southern Fe I Ka peak. The evidence for fluorescent Fe I emission from this cold gas is 
marginal at most, with Fe Ka emission extending to the north of the southern nucleus at 
3 a significance level. 



3.5. Spectral Fit with Thermal Plasma Model 

To extract the physical parameters of the highly ionized gas responsible for the Fe XXV 
and Fe XXVI emission, we fitted the 5.5-9 keV spectra of the r < 10" region. The contri- 
bution from the lower temperature gas (e.g., kT = 0.7 keV and kT = 1.4 keV components; 
Boiler et al. 2003) found in previous studies is negligible in this energy range. We adopted 
an absorbed single temperature thermal plasma {APEC; Smith et al. 2001) model for the 
collisionally ionized gas producing the highly ionized Fe XXV and Fe XXVI emission lines, 
plus zero-width Gaussian lines centered at 6.400 keV and 7.058 keV (rest frame) as the 
neutral Fe Ka and K/3 lines, respectively. A photoionized origin for the Fe XXV was not 
co nsidered as th i s was previously ruled out by the detailed photoionization calculation given 



by iNetzer et al.l (|2005[ ). The continuum contribution from the two nuclei is accounted for 
by including an absorbed power law component, representing the reflected nuclear emission. 
The absorption column {Nh = 4 x 10^'^ cm~^), the photon index (F = 0.9) and the normal- 
ization {A = 1.5 X 10~^) are fixed to be the best fit values co mbining the two nuclei. The 



abundance Z is fixed at solar values (iGrevesse fc Sauvallll998l ). 



The spectrum with the best-fit model is plotted in Figure [121 which is well represented 
with a x^/d-o.f = 300/246. The temperature of the hot plasma is 6.15 ± 0.33 keV with an 
absorption column of Nh = 5.5 ± 1.7 x 10^^ cm~^. The observed 0.5-8 keV luminosity is 
5.3 ±1.0x10^^ erg s~^ and i^o.s-skev = 5.5 ±1.4x10^^ erg s~^ after correction for absorption. 
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The normalization is 2.5 ±0.5 x 10~^, in units of [10"^'*/(47r(DA(l + z))^) / nn dV], where 
Z^A is the angular size distance to NGC 6240 (cm), rig is the electron density (cm~^), and 
uh the hydrogen density (cm~^). We obtain an emission measure (E.M.= j neUH^jdV) of 
3.1 X lO^^rj cm~^, where r] is a volume filling factor. 

Using the Fe XXV line emission as a tracer for the extent of the hottest thermal com- 
ponent, which is confined within the central r = 10", we can obtain the volume of the 
thermal gas emitting hard X-rays, V = 1.5 x 10^^?7 cm^, assuming a spherical geometry. 
This implies = 0.15ri~^ cm~^. Following the calculations done in Richings et al. (2010), 
the 70 million K gas contains a total mass of M^ot = 1.8 x IO^tj^Mq and a thermal energy 
of Eth = 6.5 X lO^^T]^ ergs. Since we assumed solar abundance, this implies the total iron 
mass in the highly ionized plasma 4.6 x IO^ti^Mq. Note that we also explored fitting with 
abundance left free and the best fit model gives Z = 0.5^qIZq with a lower kT = 4.9 ±1.1 
keV, which has minor impact on the above estimated values. 

The filling factor of the hot gas is poorly constrained observationally. Nevertheless, 
hydrodynamic simulations by Strickland & Stevens (2000) suggest that rj ranges between 
0.1% to 10% for X-ray bright galactic winds. In the starburst region t] is considered to 
be higher than in the large scale wind, close to 100% (Strickland & Heckman 2007). To 
facilitate our discussion, we assume r/ ~ 1% and caution the large uncertainty on this value. 
We consider that this is reasonable average value for the large volume spanned by the hot 
gas; the true filling factor is most probably close to unity in the innermost X-ray brightest 
region and -C 1% in the outer sphere. However, we also note that the dependence on rj is 
relatively weak for the hot gas mass and thermal energy (^72). 

4. Discussion: Origin of the Diffuse Hard X-rays 

We have unambiguously detected diffuse hard (5.5-8 keV) X-ray emission in the central 
region of the merging galaxy NGC 6240. The Fe Ka and K/3 fluorescent lines in NGC 
6240 have been firmly associated with the photoionized gas of high column density and low 
ionization close to the AGN (e.g., Komossa et al. 2003, Netzer et al. 2005). Most of the 
observed Fe I line flux is explained by the emission from the NGC 6240 nuclei, and the 
spatial concentration at the nuclear positions is consistent with this interpretation. Even 
with the long exposure, we do not find strong evidence that the dense molecular clouds in 
the nuclear region are producing fluorescent Fe lines. The excess of Fe I Ka counts in the 
CO bridge between the nuclei (Figure [TTj) is marginal {3a significance) when the scattered 
counts from both AGN are excluded. 
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What is the origin of the highly ionized gas emitting the Fe XXV hne that extends over 
5 kpc in NGC 6240? Given the LIRG nature of NGC 6240 characterised by a high SFR (25 
Mq yr~^, Engel et aL 2010), we first estimate the contribution to hard X-ray emission by the 
ensemble of X-ray binary (XRB) systems in the starburst. Using the correlation between the 
galaxy-wide 2-10 keV luminosity and the SFR provided by the Chandra survey of LIRGs 
(Lehmer et al. 2010, Table 4), the expected Lf^i^ from the XRB population is 3.3 x 10^*^ 
erg s-i in the 2-10 keV band, adopting a SFR ~ 25 yr^^ (Engel et al. 2010). This is 
only 5% of the observed 2-10 keV luminosity of the extended hard X-ray emission, strongly 
indicating that the diffuse emission is dominated by hot gas. Such hard X-ray emission is 
most likely associated with the thermal gas from merged supernova (SN) ejecta and stellar 
winds present during a starburst. The He-like Fe XXV li ne emission is also observed in other 



well k nown starburst galaxie s like NGC 253 and M 82 (e.g. lPietsch et al.ll200ll : lMitsuishi et al. 



2011 ). the ULIRG Arp 220 ( Iwasawa et al. 2005[) . and the integrated spectrum of a number 
of LIRG/ULIRG systems Jlwasawa et al.l boooh . suggesting the existence of similar high 
temperature plasma. It is worth noting that the observed diffuse hard X-ray emission in NGC 
6240 here is almost 100 times more lum inous than that observed in the classic superwind 
system M82 {L 



2-8keV 



4 X 10^^ erg s ^; 



Strickland fc Heckmaru 120071 ). 



Using optical imaging and spectroscopy to identify outflow in the warm ionized medium, 
Heckman et al. (1990) have shown NGC 6240 to host a superwind. Superwinds are believed 
to be driven by the thermal and ram pressure of an initially very hot (T ~ 10® K), high- 
pressure {P/k ~ 10'' K cm~^) and low-density wind. According to the starburst driven 
superwind model, a hot gas bubble of internally shocked wind material with a temperature 
of several keV forms in the region of intense starformation (Chevalier & Clegg 1985, Suchkov 
et al. 1994, Strickland & Heckman 2007); this hot gas eventually flows outward as a high- 
speed (few 1000 km s~^) wind. Here we detected and spatially constrained for the first time 
thermal X-ray emission from such hot (T ~ 70 MK) gas in NGC 6240, likely the thermalized 
SN ejecta within the vicinity of the starburst region. 

To further investigate whether the thermal energy content of the hot gas could be 
powered by thermalization of SNe shocks, we need to estimate the kinetic energy input 
from the SNe during the starburst. The resulting heating by SNe shocks should at least 
be comparable to Efh, considering that there are lower temperature phases and the gas has 
also bulk kinetic energy. The most often quoted SN rate for NG C 6240 is ~ 2 yr^^ (e.g.. 



van der Werf et al.lll993l : iBeswick et al.ll2001t iPoUack et al.l 120071 ) . which will be used next, 
although we note there is a small range of estimates. Engel et al. (2010) calculated a lower 
SN rate of 0.3 yr~^ assuming 20 Myr of continuous star formation at 25 Mq yr~^, based on 
the measured K-band l uminosity. Using t he co rrelation between [FeH] 1.257 /xm and the SN 
rate recently found by [Rosenberg et al.l (120121 ) from a sample of nearby starburst galaxies 



- 12 - 



(Equation 2 and Figure 9 therein), we find a higher but comparable SN rate (~3 yr^^). 
Assuming a fraction [a = 0.1, Chevalier & Clegg 1985, Thornton et al. 1998) of the kinetic 
energy input (10^^ ergs per SN) is converted into thermal energy of the hot gas, the tota l 
energy deposited during the past AtsB ~ 20 Myr of most recent starburst (ITecza et al.ll2000l ) 
is E = 4 X 10^^ ergs, which is comparable to Eth = 6.5 x 10^^ ergs. Adopting a 10% filling 
factor for the hot gas results in a larger Eth and a larger deficiency {E < Eth for a = 0.1), 
but it does not impair our conclusion, considering that the thermalization efficiency can 
be higher than the adopted v alue of 0.1 in an environment where many SNe are exploding 
(jStrickland &: Heckmanll2009l ). which will yield a larger E that is still consistent with Eth- 
These estimates demonstrate that the observed diffuse thermal gas traced by the highly 
ionized iron line emission in NGC 6240 is consistent with being heated by SNe shocks in the 
starburst. The cooling time of this thermal gas is Tcooi = Efh/Lx ~ 40 Myr, slightly larger 
than the starburst duration, which implies that such X-ray bright phase will be short lived 
once the starburst ceases. 



Following iMitsuishi et al.l ( 1201 ll ). another useful check is the comparison between the 
ejected iron mass from SNe and the o bserved iron rnass. A dopting M = 8.4 x 10~^Mq as 
the ejected iron mass per Type II SN Jlwamoto et al.lll999h and the average SN rate of ~ 2 
yr~^, the expected iron mass for the 20 Myr starburst is Mpe,sB = 3.4 x IO^Mq, which can 
easily account for the observed Mpe,x-ray = 4.6 x lO^M© considering not all the iron ejecta 
are currently X-ray emitting. Therefore this comparison also supports the above scenario 
that starburst is the origin for the 70 MK gas. 

However, the morphological comparison showed some interesting correspondences and 
differences: (1) in the nuclear region, while Fe XXV and [Fell] both peaks (Figure [9]) at 
the radio nuclei, bright H2 emission is observed between the nuclei (Figure [TOj) : (2) H2 and 
the X-ray Fe XXV emission share similar morphology at the 'arms' extending southeast and 
southwest of the southern nucleus (Figure [TU]) . which is not seen with the [Fell] emission. 
Why do the near-IR H2 and the X-ray Fe XXV emission share similar extended morphology, 
but not [Fell]? We suggest that these differences can be attributed to two aspects. First, 
there are presence of shocks in different velocity ranges. Secondly, the physical scales of SNR 
is much localized compared to starburst- driven winds. 



van der Werf et al.l ( 1l993l ) suggested that the [Fell] emission in NGC 6240 may origi- 



nate in shocks associated with SNRs during the starburst in the nuclei, based on the high 
[FeII]/Br7 ratio. When dust grains are destroyed by fast shocks, the Fe abundance in the 
gas phase is enhanced. Thus it is straightforward to expect that the [Fell] and the highly 
ionized Fe emission peaks at the nuclei of NGC 6240, where the intense star fo rmation ac- 
tivities are and young SNR are produced (ITecza et al.ll2000l : iPoUack et al.l 120071 ). Note that 
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Fe XXV in the southern nucleus is ~3 times brighter than the northern one (see Figure [TO] 
and Table [2]), in agreement with a ratio of ~2 found in [Fell] and implying a higher SN rate 
in the southern nucleus. 

From the X-ray perspective, the gas is heated to the high temperature producing the 
observed hard X-ray emission by the thermalization of the kinetic energy from an ensemble of 
SNe and winds from massive stars, inside the starburst region, and by internal shocks in the 
superwind at larger radii (Tomisaka & Ikeuchi 1988, Norman & Ikeuchi 1989, Heckman et al. 
1990). Fe XXV emission also peaks at the the radio nuclei and the peaks of [Fell] emission 
(Figure [9]). We can assume that we observe the shocked wind matter, whose temperature 
then gives us a measure of the shock velocity. For non-radiative, strong shocks in a fully 
ionized monoatomic gas, the post-shock temperature is Tgh = 3fiv'^f^/16k, where is the 
mean mass per particle and k is the Boltzmann constant (McKee & Hollenbach 1980). For 
the observed hot gas {kT = 6.15 keV, or ~ 7 x 10'' K), we obtain a shock velocity Vsh = 2200 
km s~^. Other than shocks driven by the starburst, such high velocity shocks are unlikely to 
be induced by the merger in NGC 6240: even in direct colli sion systems like th e Taffy Galaxy 
(UGC 12914/5), the progenitors collided at ~600 km s~^ (IBraine et al.ll2003l ). The merging 
galaxies in NGC 6240 have well passed the first encounter, without showing morphology 
like the ring galaxies exemplifying such a strong direct collision, and the estimated collision 
velocity from the CO line width is only between 150 and 300 km s~^ (IWang et al.l Il991 
Tecza et al.lboooh . 



In contrast, the bright H2 emission observed between the nuclei of NGC 6240 has been 



interpreted to originate in gl obal dissipative shocks resulting from the co l 



of the merging disk galaxies (IHerbst et al.l Il990l : Ivan der Werf et al.l 



The shock velocity was constrained to be at most 40 km s ^ (Ivan der Werf et al 



1993 



isions of the ISM 



Max et al. 



2005). 



19931). 



However, the fainter large scale 'arms' of H2 emission extending SE and SW of the southern 
nucleus needs a different explanation, as it closely follows the X-ray emission. As fast 
outflows (several thousand km s~^) driven by the starburst are expected (e.g., Heckman et 
al. 1990, Fabbiano et al. 1990; see Veilleux et al. 2005 for a review), this large scale H2 
emission can be interpreted as originated from molecular material entrained and shocked by 
the superwind. The thermal velocity of this 6 keV gas, estimated as the adiabatic sound 
speed Vth = (5fcT/3/xmH)^''^ ~1260 km s~^ (assuming solar abundance), is much higher than 
the escape velocity Vesc = \/ '^GM[t, which is ~2 10 km s~^ at r ~ 1 kpc using the dynamical 
mass Miiyn = 6 x lO^M© in iTacconi et al.l ( 1l999l ). Therefore the expanding hot bubble can 
"blow out" to larger radial distance such as the 5 kpc exten t seen here. Following a s imple 
model of a shock propagating in an inhomogeneous ISM in Ivan der Werf et al.l ( 1l993l ). the 
high velocity shock {yx = 2200 km s~^, see above) in the low density wind responsible for the 
Fe XXV emission will slow down significantly when encountering a layer of cold molecular 
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gas. The decreased shock speed V2 can be estimated using the density contrast between the 



X-ray gas (pi ~ O.lr/"^/^ cm~^) and the molecular cloud (p2 ~ 10^ ~ 10^ cm~^), following 
^2/^1 = a/pi/p2 (Spitzer 1978, Equation [12-39]; note this is different from the shock jump 
condition). The expected f 2 ~ 7 — 20 km s~^ agree s well with the velocity c onstraint from 



shock excitation for the observed H2 emission (e.g., Ivan der Werf et al.l 119931 ). 



More evidence supporting this scenario also comes from the presence of high velocity H2 
and Ha emission in the 'arms' region, based on prominent wings in the H2 1-0 S(l) line (with 
full width at ze ro intensity of ^1600 k m s~^; van der Werf et al. 1993; Engel et al. 2010) and 
Ha line profile ( iHeckman et al.lll990[ ). tracing the ambient gas entrained and shocked by the 
superw ind. In a very recent CO(l-O) interferometry observation of NGC 6240. lFeruglio et al. 
(120131 ) further identified evidence for a shock wave that propagates eastward from the nuclei. 
One of the resolved structures of blue shifted CO emission is spatially coincident with the 
SE H2 filament (Max et al. 2005), which is also associated with Ha and soft X-ray emission. 
The derived kinetic power of the outflowing CO gas is 6 x lO''^ erg s~^, with an estimated 
age > 2 X 10'' years. This implies a total kinetic energy of = 3.8 x 10^'^ ergs. Recall that 
the available kinetic input from the multiple SNe explosions during the starburst is 4 x 10^^ 
ergs (10^^ ergs per SN x 2 yr~^ x2 x 10^ yr), there seems no deficit in the energy budget. 
However, we cannot yet draw conclusion that no additional energy injection is required, 
until a full tally of the bulk kinetic energy and thermal energy of the outflowing gas in other 
phases is available for NGC 6240. Under the assumption that the 70 million K gas is volume 
filling [r] = 100%; Strickland & Stevens 2000), its Eth would exceed the available kinetic 
energy provided by the recent starburst by 2.5 x 10^^ ergs, implying additional kinetic power 
input such as from an AGN outflow or multiple star forming episodes. In fact, the thermal 
energy cont ent in the larg e scale (~100 kpc) soft X-ray halo {kT = 0.65 keV) in NGC 6240 
identified b yNardini et al.l (l2013[ l is estimated to be 4.9 x 10^*^ ergs, which hints at additional 
energy inp ut from a widesprea d, enhanced star formation proceeding at steady rate over 
~200 Myr ( iNardini et al.ll2013l ) besides the most recent nuclear starburst. 



Lastly, we note that a good spatial correlation is al so expected if the H2 emission is 
excited by localized X-ray irr adiation of molecular clou ds ( iDraine fc Woodslll990l ). However, 
this was concluded unli kely ( Ivan der Werf et al.lll993l ). since the required SN rate of 7 yr~^ 
( IDraine fc Woodslll99ll ) is much higher than the inferred rate from radio, infrared, and our 
X-ray measurement. Thus we conclude that propagation of fast shocks originated in the 
starburst driven wind into the ambient dense gas can account for the large scale X-ray and 
H2 morphological correspondence. 
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5. Conclusions 

Combining new and archival Chandra observations of NGC 6240, we have obtained the 
deepest X-ray image of the central r = 5 kpc (10") of the merging ULIRG NGC 6240 with 
sub-arcsecond resolution. Studying the hard X-ray emission centered at Fe K complex and 
comparing to images in other wavebands, here for the first time we clearly resolved hard 
extended emission with its extension established, and assessed the nuclei contribution to the 
hard emission separately. Such a hard extended emission was predicted to be the innermost 
manifestation of superwinds. So far, the superwind nature of NGC 6240 had been studied 
mostly in the Ha, molecular gas, and soft X-rays. We summarize our findings as follows: 

1. The X-ray emission from the pair of active nuclei and the diffuse hot gas is resolved 
in detail. Extended X-ray emission from kT ~ 6 keV (~70 million K) hot gas was 
found over a spatial scale of 5 kpc, indicating the presence of fast shocks with velocity 
of ~2200 km s-^ 

2. The observed luminosity of the 70 million K temperature gas is i^o.5-8kcV = 5.3 x 10^^ 
erg s~^, with a total mass of M^ot = 1-8 x lO^M© and thermal energy Et^ = 6.5 x 10^''' 
ergs, assuming a filling factor of 1% for the hot gas. The total iron mass in the highly 
ionized plasma is Mpexxv = 4.6 x lO^M©. Both the energetics and the iron mass in the 
hot gas are commensurate with the expected injection from the starburst, assuming 
an average SN rate of 2 yr~^ that is in agreement with its high star formation rate and 
estimates from non-thermal radio emission and [Fell] luminosity. 

3. The spatial distribution of the highly ionized gas emitting Fe XXV peaks together with 
the [Fell] emission around the radio nuclei, associated with the high concentration of 
young SNRs during the nuclear starburst. The extended Fe XXV emission shows a 
remarkable correspondence to the large scale morphology of H2(l-0) S(l) line emission 
and Ha filaments. Propagation of fast shocks originated in the starburst driven wind 
into the ambient dense gas can account for this morphological correspondence. The Fe 
I Ka emission peaks at the active nuclei. No evidence for fiuorescent Fe I emission is 
found in the CO filament connecting the two nuclei. 

Although the above superwind scenario accounts for the extended hard X-ray emission 
and observational evidence in other wavebands, we cannot yet rule out additional energy 
injection in the central region of NGC 6240 (e.g., heating by an AGN outflow). Given the 
poorly constrained volume filling factor of the X-ray emitting gas, the thermal energy of 
the hot gas may exceed the available kinetic energy provided by the SNe explosions during 
the most recent starburst (Efh — Esb — '^■^ x 10^^ ergs for 77 = 100%), implying additional 
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energy input such as from an AGN outflow or more complex star formation activities besides 
the recent nuclear starburst (see discussion in Nardini et al. 2013). A full tally of the bulk 
kinetic energy and thermal energy of the outflowing gas in other phases is crucial for the 
evaluation of energy budget in the NGC 6240 system. We will further characterize the soft 
X-ray emission within the central 15 kpc in this merging galaxy in a forthcoming paper 
(Wang et al. in preparation). 
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Table 1: Chandra ACIS Observation Log for NGC 6240. 



ObsID 


Date 


Exposure (ks) 


Instrument 


Data Mode 


Roll Angle (°) 


1590 


Feb-22-2002 


37 


ACIS-S 


Faint 


246.2 


6909 


May- 11-2006 


143 


HETG/ACIS-S 


VFaint 


129.9 


6908 


May-16-2006 


159 


HETG/ ACIS-S 


VFaint 


138.1 


12713 


May-31-2011 


146 


ACIS-S 


Faint 


167.2 


Total Effective Exposure* (ks) ... 




363 









* The effective exposure time for the Chandra HETG zeroth order image refers to the equivalent exposure 
needed if the data were taken with direct Chandra ACIS imaging, calculated by comparing its effective area 
to that of ACIS imaging without grating. 
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Table 2: Chandra ACIS Measurements of the Complex Fe K Lines in NGC 6240. 



Region 


Counts (5.5-8 keV) 




Fe I 


Fe XXV 


Fe XXVI 






Energy (keV) * 


6.26±0.01 


6.54±0.03 


6.86±0.04 


Nuc-N 


908 


Flux ** 


10.7±0.8 


2.2±0.5 


1.6±0.6 






EW (eV) 


645158 


731:11 


105+32 






Energy (keV) * 


6.23±0.01 


6.46±0.01 


6.681-- 


Nuc-S 


2096 


Flux ** 


13.5±0.9 


8.0±0.8 


4.5±0.7 






EW (eV) 


342±30 










Energy (keV) * 


6.25±0.01 


"■^^-0.02 


6.74±0.02 


r < 10" 


4856 


Flux ** 


28Atli 


16.2±1.1 


4.3t?-i 






EW (eV) 


359±22 




53±14 






Energy (keV) * 


6.26±0.02 


6.52±0.02 


6.84±0.04 


r < 70" 




Flux ** 


26±8 


18±7 


9±7 


(XMA/)t 




EW (eV) 


300±100 


220±90 


120±90 



* Observed energies of the line centers. 

** Photon flux in the unit of 10~^ photons s~^ cm"^. 

XMM- Newton EPIC measured FeK line complex from Boiler et al.(2003). Because of the lower angular 
resolution of XMM-Newton, the spectrum is extracted using a circular region with radius of 70". 
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Fig. 1.— (a) HST/ACS F814W image of the central 25" x 25" (12 kpc across) region of NGC 
6240, showing the disturbed optical disk, (b) The Chandra ACIS image (5.5-8 keV) of the 
same region extracted from the merged data. The white circles (r = 10") mark the spatial 
extent of the hard X-ray emission. The green ellipses (4" x 2") enclose the two nuclei. 
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Fig. 2. — Comparison between the observed surface brightness profile (black) and that 
expected from the Chandra PSF for the 5.5-8 keV emission (red). Note the presence of 
significant hard X-ray emission reaching r ~ 10" radial distance. 
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Fig. 3. — Comparison between the X-ray spectra for the dual AGN region (shown in blue; 
extracted from the ellipse in Figure [1]) and the 10 arcsec radius region (shown in black; 
extracted from the circle in Figure [1]) including both the diffuse and AGN emission. The 
X-ray spectrum for the central region excluding the nuclei is also shown (green) to illustrate 
the spectral difference and its flux level. The energy scale is observed (not rest frame). 
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Energy (keV) 




Energy (keV) 




Fig. 4. — Top to bottom: ACIS spectra and fits for tfie extended central emission, tfie 
nortliern nucleus, and the southern nucleus. Both ACIS data (black) and HETG (red) 
zeroth order data were used. The model components are also shown as dashed lines (see 
Table E]). 
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Fig. 5. — Narrow band ACIS images: (a) the 5.5-6 keV continuum band; (b) Fe Ka line (6-6.4 
keV); and (c) Fe XXV line (6.4-6.7 keV). The latter two images were continuum-subtracted 
(see text). 




Fig. 6. — Narrow band ACIS images, with EMC2 PSF-deconvolution applied: (a) the con- 
tinuum band (5.5-6 keV); (b) hne image for the neutral Fe Ka line (6-6.4 keV); and (c) the 
FeXXV (6.4-6.7 keV). 
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Fig. 7. — 2-D image fitting of tlie innermost nuclear region in (a) 5.5-8 keV (b) 5.5-6 keV 
(continuum) (c) 6.0-6.4 keV (Fe I) (d) 6.4-6.7 keV (Fe XXV). For eacli band, tlie observed im- 
age, simulated pair of point sources representing tfie double nuclei, an d the residual emiss ion 
are shown. The two circles mark the compact radio nuclear sources (IColbert et al.lll994j ). 
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Fig. 8. — (a) The adaptively smoothe d 5.5-8 keV ACTS image (in 1/2 pixel binning) overlaid 
with the HST/WFPC2 Ha contours Jhira et al.l[2nn3 ). (b) The PSF-deconvolved 5.5-8 keV 
image overlaid with the same Ha contours, (c) The PSF-deconvolved continuum band (5.5- 
6 keV) image (Figure E^) overlaid with the same Ha contours, (d) The PSF-deconvolved 
Fe XXV image (Figure ^) overlaid with the same Ha contours. 




Fig. 9. — The near-IR [Fell] 1.64 /iin contours (van der Werf et al. 1993) overlaid on: (a) 
raw (6-6.4 keV) Fe I Ka emission; (b) raw (6.4-6.7 keV) Fe XXV line emission; (c) same as 
(a) but with EMC2 deconvolution; (d) same as (b) but with EMC2 deconvolution. 
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Fig. 10. — (a) The FeXXV line emission image overlaid with the contours of H2(l-0) 2.12 /im 
emission (Max 2005); (b) PSF-deconvolved Fe XXV line emission image with the same H2 
contours. 
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Fig. 11. — (a) The Fe Ka line emission image overlaid with the contours of SMA CO(3-2) 
emission (U et al. 2011); (b) PSF-deconvolved Fe Ka line emission image with the same CO 
contours. 
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Fig. 12. — The X-ray spectra of the NGC 6240 r = 10" central emission (black: ACIS data; 
red: HETG zeroth order data) and the fit using thermal plasma model for the hot gas (see 
the text in § 13. 5p . 



